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The  meso-scale  flow  over  and  around  a  finite  lengtn  ^ moun'Ain*ridge  iaftdtjl  U 
typical  mountain  valley  terrain  is  studied.  The  present  study  is  different 
from  earlier  analytical  lee  wave  studies  in  that  the  solution  is  truly  three- 
dimensional  for  more  realistic  topography.  A  simple  model  is  developed  and  is 
attacked  by  multiple-scaling  approach.  The  solutions  are  not  only  valid  for 
large  distances  in  the  downwind  direction,  but  are  also  good  for  the  near 

regions.  The  phase  lines  of  the  vertical  velocity  on  a  horizontal  plane  at  a 
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given  height  have  shape  of  hyperbolas  and  are  concave  toward  downwind  side  in 
accord  with  observations.  The  wind  component  parallel  to  the  mountain  ridge  is 
found  to  be  diffluent  on  the  windward  slope  side,  continuing  to  be  so  after 
crossing  over  the  ridge  in  the  lowest  layer,  and  to  become  confluent  at  some 
distance  downwind  from  the  ridge.  The  intensity  of  the  lee  wave  over  a  valley 
relates  to  the  separation  between  the  two  ridges  and  the  stability.  Just 
above  the  valley  floor,  v/ind  parallel  to  valley  axis  is  in  the  down  valley 
direction  and  thus  raises  the  height  of  the  maximum  down-valley  wind  much 
higher  than  that  of  maximum  'slope  wind'  above  the  valley  sides  during  the  j 
night.  The  computation  also  shows  the  increase  of  its  intensity  and  height  of 
occurance  with  down-valley  distance.  In  the  farther  downwind  direction  of  the 
lee,  the  maximum  confluent  flow  along  the  ridge  is  found  at  a  much  lower  level 
and  is  stronger.  The  intensity  and  the  position  of  the  maximum  viind  parallel 
to  the  valley  or  the  mountain  ridges  is  closely  related  to  the  downslope 
motion  of  the  lee  wave. 
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prefacf: 


The  thermall  y  imluced  mean  t  a  in -valley  '/ind  system  and  its 
interaction  with  the  large  or  lua-uo-scale  flow  Jms  always  inter¬ 
ested  meteorologists  and  other  professionals  directly  involved 
with  the  development  of  the  wind  svstem.  VarioiiS  analytical  and 
numerical  studies  have  been  pursued  and  have  yielded  some  inter¬ 
esting  results.  Yet  many  phenomena  in  the  local  wind  system  liave 
not  really  been  explored  and  explained  satisfactorily  F'oi' 
example,  the  problem  concerning  the  speed  and  lieight  of  tlie  maxi¬ 
mum  wind  or  local  jet  along  the  \/dliey  axis  I'emains  unclear. 

Whether  they  result  from  direct  influence  oi  the  basic  background 
flow  or  just  the  f)ure  local  theimial  boundary  larger  effect  or  from 
interaction  of  both  requires  a  further  understanding  of  some  basic 
problems  of  flow  over  and  around  typical  mountain-valley  terrain 
and  the  development  of  a  better  method  to  solve  the  complicated 
problem  of  sophisticated  model  of  boundary  layer  dynamics. 

In  the  first  part,  we  developed  a  model  to  study  the  flow 
over  and  around  more  realistic  three-dimensional  mountain-valley 
terrain  and  utilized  a  multiple  scaling  approach  to  obtain 
truly  three-dimensional  flow  solutions  for  the  near  region,  as  well 
as  for  a  large  distance  from  the  ridges  in  the  lee.  He  theoretical 
findings  from  the  lee-wave  over  the  lee-&lojj>e  are  in  accord  with 
some  observations  and  other  studies.  This  investigation  is  dif¬ 
ferent  from  many  other  previous  studies  in  that  the  flow  around 
the  ridge  and  along  the  valley  has  enabled  us  to  visualize  how 
interaction  between  the  meso-scale  flow  and  the  valley  wind  system 
in  the  lower  layer  takes  place  and  to  understand  better  the  wind 
structure  u  .der  such  circumstances.  The  solution  for  the  wind 
component  xn  a  direction  parallel  to  the  ridges  or  in  the  valley 
direction  explains  some  puzzling  and  interesting  valley-wind 
phenomena . 
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For  the  second  part,  studies  of  boundary  layer  problems 
are  presented.  The  parameterization  of  the  surface  layer  and 
active  soil  layer  with  moisture  and  precipitation  contributions 

I 

is  made  for  meso-scale  dynamics .  It  can  be  used  for  mountain- 
valley  circulation  as  long  as  the  slope  of  the  terrain  is  not 
exceedingly  large  and  some  other  basic  conditions  are  met. 

With  this  parameterization  the  boundary  conditions  for  the 
governing  prediction  equations  of  mountain-valley  wind  problems 
in  numerical  models  can  be  appropriately  constructed.  Applying 
these  improved  boundary  conditions  will  facilitate  the  solution 
of  the  numerical  model  of  the  boundary  layer  and  over  meso-  or 
large-scale  circulation  problems  and  can  also  produce  more  real¬ 
istic  results. 
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ABSTRACT 


The  meso-scale  flow  over  and  around  a  finite  length  of  mountain  ridge 
and  typical  mountain  valley  terrain  is  studied.  The  present  study  is  dif¬ 
ferent  from  earlier  analytical  lee  wave  studies  in  that  the  solution  is  truly 
three-dimensional  for  more  realistic  topography.  A  simple  model  is  developed 
and  is  attacked  by  multiple-scaling  approach.  The  solutions  are  not  only 
valid  for  large  distances  in  the  downwind  direction,  but  are  also  good  for  the 
near  regions.  The  phase  lines  of  the  vertical  velocity  on  a  horizontal  plane 
at  a  given  height  have  shape  of  hyperbolas  and  are  concave  toward  downwind 
side  in  accord  with  observations.  The  wind  component  parallel  to  the  mountain 
ridge  is  found  to  be  diffluent  on  the  windward  slope  side,  continuing  to  be  so 
after  crossing  over  the  ridge  in  the  lowest  layer,  and  to  become  confluent  at 
some  distance  downwind  from  the  ridge.  The  intensity  of  the  lee  wave  over  a 
valley  relates  to  the  separation  between  the  two  ridges  and  the  stability. 

Just  above  the  valley  floor,  wind  parallel  to  valley  axis  is  in  the  down  valley 
direction  and  thus  raises  the  height  of  the  maximum  down-valley  wind  much  higher 
than  that  of  maximum  'slope  wind'  above  the  valley  sides  during  the  night.  The 
computation  also  shows  the  increase  of  its  Intensity  and  height  of  occurance 
with  down-valley  distance.  In  the  farther  downwind  direction  of  the  lee,  the 
maximum  confluent  flow  along  the  ridge  is  found  at  a  much  lower  level  and  is 
stronger.  The  intensity  and  the  position  of  the  maximum  wind  parallel  to  the 
valley  or  the  mountain  ridges  is  closely  related  to  the  downslope  motion  of  the 
lee  wave. 
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I ,  INTRODUCTION 


Previous  theoretical  studies  of  ]ee  waves  are  numerous.  Most  attention 
has  been  focussed  on  two-dimensional  lee  wave-waves  behind  an  infinitely  long 
ridge.  Earlier  works  on  three-dimensional  lee  waves  are  essentially  limited 
to  simple  topography. 

Scorer  and  Wilkinson  (1956)  studied  waves  in  the  lee  of  both  Isolated 
circular  and  oval-shaped  hills  in  a  two-layer  model  atmosphere.  In  each  layer, 
the  "Brunt-VMisMia  frequency  over  mean  wind  speed"  is  assumed  constant.  They 
found  that  the  wave  disturbance  is  essentially  confined  to  a  wedge  (for  example, 
with  wake  angle  of  24°  in  one  case)  and  has  the  form  of  a  shipwave.  Assuming 
an  isothermal  atmosphere  with  uniform  wind  velocity,  and  ground  surface  as  a 
stepped  plateau  of  finite  width,  Wurtele  (1957)  was  able  to  show  the  wave 
nodal  lines  in  a  horizontal  plane  are  rectangular  hyperbolas  concave  downwind, 
in  qualitative  agreement  with  clouds  observed  in  the  lee  of  Mt.  Fuji  and  lab¬ 
oratory  experiments  by  Abe  (1941).  Although  he  cannot  use  his  one  layer  model 
to  show  the  shipwave  type  disturbance,  Wurtele  pointed  out,  in  a  later  paper 
concerlng  cloud  photograpnies  for  the  lee  of  San  Nicolas  Island  (Endinger  and 
Wurtele,  1972)  that  when  a  strong  temperature  inversion  occurs,  a  classical 
shipwave  disturbance  will  develop  instead  of  a  cresent-shaped  pattern. 

Palm  (1958)  considered  the  case  when  the  basic  wind  increases  linearly 
with  height  in  the  troposhere  and  stays  constant  in  the  stratosphere,  while  the 
stabilities  are  constant  in  both  layers.  For  a  circular  hill.  Palm  obtained  a 
solution  for  the  vertical  velocity  and  found  that  the  lee  waves  are  also  confined 
to  a  wedge  in  the  downwind  direction,  but  only  for  sections  at  a  relatively  large 
distance  from  the  hill.  In  extending  his  method  of  study  of  two-dimensional 
lee  waves,  Onlshl  (1960)  sought  a  solution  for  three-dimensional  lee  waves  due 
to  a  circular  hill  in  the  atmosphere  with  both  zero  and  constant  wind  shear  in 
the  vertical  direction.  As  done  by  Scorer  and  Wilkinson,  Wurtele,  and  Palm, 
Onlshl  also  applied  the  method  of  stationary  phase  to  obtain  the  solution.  For 
the  case  of  constant  shear,  he  managed  to  increase  the  wake  angle  from  a  small 
wedge  to  a  value  beyond  which  no  stationary  point  can  be  obtained. 
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In  consequence,  there  Is  no  iee  wave  found  outside  this  maxitium  wake 
angle.  As  for  the  case  where  there  is  constant  wind  and  strong  stability, 
the  disturbance  is  confined  to  a  strip.  Crappe?  investigated  lee  wave 
problems  for  a  circular  mountain  first  and  then  a  mountain  with  elliptical 
contours  (1962).  For  the  latter  study,  constant  wind  shear  together  with 
both  thermal  stability  conditions  constant  and  an  exponentially  decreasing 
with  height  were  considered.  Under  the  constant  stability  condition,  a- 
found  that  the  waves  are  confined  to  a  strip  for  the  case  of  no  shear, 
whereas  the  waves  are  confined  to  a  wedge  resembling  shlpwaves  for  the  case 
with  a  large  constant  shear.  For  the  case  without  wind  shear,  the  amplitude 
of  the  waves  in  the  central  plane  in  the  downwind  direction  of  an  elliptical 
mountain  is  greater  than  that  of  an  infinite  ridge.  The  magnitude  of  the 
vertical  disturbance  decays  in  a  form  inversely  proportional  to  three  halves 
(3/2)  power  of  the  downwind  distance.  The  results  appear  to  be  controversial. 

Evaluation  of  integrals  by  numerical  methods  or  the  use  of  numerical 
models  for  studies  of  three-dimensional  lee  wave  problems  has  Increased 
recently.  A  three-dimensional  lee  wave  for  an  isolated  nearly  circular  hill 
was  studied  with  the  perturbation  method  but  integrated  numerically  in  both 
two-  and  three-layer  models  by  Sawyer  (1962).  In  his  model,  the  horizontal 
wind  velocity  can  change  with  height  in  direction,  as  well  as  in  speed.  The 
computed  vertical  velocities  due  to  the  lee  waves  have  a  series  of  bow-shaped 
patterns  concave  downwind  and  which  are  also  confined  to  a  wedge  similar  to 
those  obtained  by  others  The  maximum  magnitude  occurs  just  off  the  peak  of 
the  hill  and  decreases  gradually  away  from  the  hill.  By  being  aware  that  the 
lee  wave  disturbance  is  essentially  confined  to  a  wedge  and  applying  a  trans¬ 
formation  to  an  oblique  co-ordinate  system,  Pekells  (1966)  developed  an 
algorithm  for  solving  the  non-linear  three-dimensional  lee  wave  problem.  Com¬ 
putations  show  that  the  height  of  an  obstacle  definitely  influences  the  shape 
of  the  pattern  of  vertical  velocity  in  the  non-linear  model.  For  high  steep 
mountains  where  the  non-linear  process  is  important,  the  pattern  differs 
strikingly  from  that  of  relatively  low  obstacles  in  that  the  nodal  lines  of 
the  vertical  velocity  are  convex  downwind,  almost  opposite  to  the  shipwave  or 
bow-shape  wave. 
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Vergelner  (1975,  1976)  cltivelop^cl  a  three-dimensional,  time  dependent, 
linearized  numberlcal  model  to  study  Poehn  wind  and  iee  wave  flow.  For  the 
circular  hill  case  he  produced  a  more  or  less  similar  kind  of  vertical  velocity 
pattern  as  obtained  by  ethers  except  that  some  computing  noises  are  Involved. 

He  also  applied  the  model  to  study  the  flow  over  the  smoothed  topography  of 
a  part  of  the  Tyrol  Province  of  Austria  and  selected  soundings  yielded 
interesting  results.  Some  topographically  produced  features  under  the  inter¬ 
action  between  orographlcally  forced  and  free  waves  are  still  recognizable  in 
the  computed  wind  field.  The  computed  results  are  generally,  also,  in  agree¬ 
ment  with  other  theoretical  and  observational  evidence. 

Considering  both  cross-wind  and  vertical  shear,  Blumen  et  al.,(1975) 
attacked  the  three-dimensional  lee  wave  behind  a  nearly  circular  mountain. 

The  vertical  velocity  field  is  obtained  in  Fourier  integeral  and  evaluated 
numerically.  One  of  the  interesting  results  is  that  when  the  constant  basis 
wind  velocity  is  23  msec  \  the  maximum  upward  vertical  velocity  is  about 
2.5  cm  sec  ^  at  height  2  km  and  distance  2  km  downwind  from  the  center  of  the 
mountain  which  has  a  peak  height  of  about  570  m.  This  vertical  velocity  is  a 
much  lower  value  than  the  one  usually  obtained  by  other  authors. 

Most  recently,  Gjevik  et  al  (1978)  applied  a  kinematic  method  developed 
by  Whltham  to  investigate  the  wake  angle  of  shipwave  type  patterns  of  the  phase 
lines  of  the  vertical  velocity.  The  computed  wave  patterns  under  various 
conditions  are  qualitatively  in  agreement  with  some  satellite  cloud  photographs. 

The  purpose  of  the  present  investigation  is  to  study  the  three-dimensional 
lee  wave  above  the  flow  over  and  around  a  ridge  of  finite  length  and  above  a 
typical  mountain-valley  terrain,  and  its  effects  on  the  mountain-valley  wind 
system.  The  problem  is  studied  based  on  a  simple  linearized  model  and  is  solved 
with  the  multiple  scaling  approach  so  that  representative  solutions  can  be 
obtained  for  both  the  near  region  and  the  far  region  of  the  mountain  terrain  in 
the  downwind  direction.  From  the  theoretical  study,  the  basic  structure  of  the 
velocity  field  in  the  lee  of  the  finite  ridge  and  above  a  typical  valley  is 
obtained  when  the  thermal  boundary  layer  influence  is  minimum.  Because  of  our 


3 


interest  in  effects  on  valley  winds,  we  direct  more  attention  to  the  near 
field  than  has  often  been  the  case  in  the  lee  wave  studies.  In  Section  2, 
a  general  description  is  presented.  The  development  of  the  solutions  based 
on  the  multiple-scaling  method  and  the  numerical  results  of  the  solutions 
are  given  in  Sections  3  and  4.  A  brief  conclusion  and  recommendations  can 
be  found  in  Section  5. 
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II.  THE  MODEL 


Since  our  main  interest  is  to  explore  the  consequences  of  mountain 
waves  for  valley  winds,  we  simplify  as  -iiuch  as  possible  the  mountain  wave 
aspect  by  considering  a  time-independent  Bousslnesq,  incompressible,  stratified 
invlscld  fluid,  and  linearize  about  a  basic  state  with  velocity  U(Z)>0  in  the 
x-dlrection,  and  thermal  stratification  T^(Z).  Thus,  our  basic  equations  for 
the  perturbation  velocity  u  and  temperature  6  are: 

Uu  +  wu'i  +  Vp  -  ag0K  =  0,  V*u  =  0,  and  US  +  wT'  *  0. 

—X  —  P  xo 

o 


(u  is  the  thermal  expansion  coefficient.)  Eliminating  all  the  dependent 
variables  except  w  amont  these  equations  one  obtains 


2 

V  w 

XX 


r'w 


XX 


~2 

U  ^ 


=  0 


(1) 


where  is  the  horizontal  Laplace  operator  and  '  is  the  square  of 

the  Brunt-vaisMia  frequency. 


We  take  this  equation  on  O^Z^irH  with  the  boundary  conditions 
w  *  U(o)  ^  on  Z  =  0  (2a) 

aX 

w  “  0  on  Z  =  ttH  (2b) 

The  function  h  in  (2a)  represents,  on  the  linear  theory,  the  height  of  the 
mountain  as  a  function  of  the  horizontal  coordinates.  We  are  going  to  assume 
that  the  characteristic  scale  of  variation  of  h,  in  its  dependence  on  x,  is  of 
the  same  order  as  the  depth  ttH  of  the  layer.  However,  its  characteristic  scale 
of  variation  in  y  is  taken  to  be  larger  by  a  factor  of  e  where  e  is  a  small 
parameter.  It  is  this  aspect  ratio  parameter  e  which  we  shall  exploit  in 
solving  the  problem;  the  small  parameter  implicitly  already  used  in  making  the 
linearization  may  be  taken  to  be  the  ratio  h^/H,  where  h^  is,  say,  the  maximum 
height  of  the  mountain. 


We  assume  h  tends  to  zero  (with  sufficient  rapidity  for  mathematical  needs) 
as  X  or  We  shall  give  explicit  calculations  for  the  example: 

h  =  ho  exp  k|x/h|  -  (ey/H)^!  (3) 

(In  place  of  (2b)  one  could  also  consider  an  unbounded  layer,  H  then  being 
an  appropriate  scale  height  for  the  velocity  and  stratification  profiles; 
if  the  latter  are  such  as  to  admit  wave  propagation  at  large  Z,  the  condition 
w->o  would  then  be  replaced  by  an  outgoing  wave  condition.) 

We  introduce  dimensionless  position  coordinates  temporarily  distinguished 
by  overbars,  by 

X  =  xH,  y  =  yH,  Z  =  zH 

and  also  set  ey  =  Y.  We  set 

h  =  hoh(x,Y),  U  =  U(o)U(z),  N^U~“  =  H*^J(i) 

and  for  the  dependent  variables  take: 

u  =  U(o)h^H  ^u,  V  =  eU(o)h^H  ^v, 

w  =  U(o)h^H“^“^  e  =  U(o)^h^(agH^)"^0 

P  =  p^U(o)^h^H  V 

Then  the  basic  equations  become,  on  dropping  the  overbars; 


Uu  +  U'w  +  p  = 

X  X 

0 

(4a) 

Vv  +  Py 

X  I 

0 

(4b) 

Uw  -6  +  P.,  = 

X  z 

0 

(4c) 

e  +  ujw 

X 

0 

(4d) 

X  2  . 

U  +  EV„+W” 

X  Y  Z 

0 

(4e) 

6 


(5) 

(6a) 

(6b) 


We  have  used  U(o)  as  the  characteristic  scale  for  the  velocity  profile; 
this  would  of  course  not  be  appropriate  if  U(o)  =  0;  by  in  such  a  case 
neither  would  be  the  linearization.  In  the  calculations  using  (3),  which 
in  dimensionless  form  is 


we  shall  also  take  U  and  J  to  be  constants. 


In  addition  to  these  boundary  conditions,  we  have  the  conditions  of 

no  upstream  influence,  meaning  that  w  (and  the  other  perturbation  quantities 

as  well)  should  tend  to  zero  for  x-*-  and  downstream  boundedness,  meaning 

that  they  should  bounded  for  x-+  +«.  We  also  assume  that  they  tend  to  zero 

as  Y->±®  for  any  fixed  x.  Actually,  we  are  considering  a  family  of  problems 

parameterized  by  e,  and  we  shall  require  these  conditions  expressed  in  the 

above  dimensionless  form,  to  hold  uniformly  in  e,  as  We  shall  also  be 

concerned  to  find  a  representation  of  the  solutions  which  is  uniformly  valid 

In  as  large  a  domain  as  possible  as  e-Kj.  If  the  problem  is  attacked  as  a 

perturbation  problem  in  e  in  the  most  straightforward  way,  one  does  not,  in 

fact,  get  a  representation  uniformly  valid  for  large  x;  in  order  to  obtain  a 

relatively  simple  and  useful  description  of  the  flow  a  somewhat  different 

approach  is  needed.  We  shal’  see  that  there  are,  in  fact,  three  different 

-1  -2 

scales  in  x,  X  “  0(1),  0(e  ),  and  0(e  ),  which  are  important  in  this  descrip¬ 

tion;  this  will  be  clarified  by  consideration  of  the  specific  example  (7), 
to  which  we  now  proceed. 
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Ill  THE  SOLUTIONS 


We  now  take  U  and  J  to  be  constants,  so  U  =  1;  we  assume  J>0  (static 
(n) 


stability).  Let  w 
0<Z<iT : 


(x,Y)  to  be  the  Fourier  sine  coefficients  of  w  on 


.(n) 


sin  nz  w  dz 


We  shall  represent  w  by  its  sine  series,  but  since  w  does  not  vanish  on 
Z  =  0,  the  terms  of  this  series  will  only  tend  to  zero  like  1/n.  Thus  the 
series  will  not  be  directly  useful  for  computation;  we  shall  consider  later 
how  it  can  be  improved.  More  important  is  the  fact  that  formal  term  by  term 
differentiation  of  the  series  does  not  give,  for  example,  the  correct  sine 
series  of  ^  •  Thus,  the  sine  series  should  not  simply  be  substituted 


into  eq  (5).  However,  the  consequences  of  (5)  for  the  w^"^  are  readily 

obtained  by  multiplying  it  by  (2/tt)  sin  nz  and  integrating  from  0  to  it.  On 

2  2 

integrating  the  term  involving  3  w/3z  twice  by  parts,  using  equations  (6) 
one  obtains: 


1  \ 

o  / 

32  2  ' 

—7  -  n  J 

— 2“  ^  ^ 

1  3x 

8x 

3 


2 

) 

9x 


9y 


Ir  Ui 

IT  .  3 
3x 


(3) 


Now  take  the  Fourier  transform  in  Y  of  this  equation,  setting  w^”^ 

—00 


this  gives: 

/  2 

3  2 


-  n  +  J| 


\3x 


3x^ 


e^k^ 


3 

— 2 
\  3x 


+  J  I  w<"> 


^  3_h 

"  3x^ 


(9) 


For  h  given  by  (7),  we  have 


tr^  I — -K|x|-k  /4  .  X 

3  h  ■  -K  vrre  sgn(x), 


3x 


(10) 
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f 


Because  of  the  discontinuity  in  ^  at  x  =  0  we  shall  consider  separately 
the  cases  of  positive  and  negative  x,  applying  appropriate  matching  conditions 
at  the  origin.  To  obtain  these  conditions,  we  look  also  at  the  other  com¬ 
ponents  of  the  flow,  representing  u,  v  and  p  by  cosine  series  of  the  form 


00 


r  (n)  (n)  2 

}  u  cos  nz,  u  =  —  \  i 

I  I  Jo 


u  cos  nz  dz I  ,  and  6  by  a  sine  series 


as  in  the  case  of  w.  Taking  also  the  Fourier  transform  with  respect  to  Y, 
equations  (4)  (for  the  present  example)  then  give: 


=  0 

X  x 


(n>o) 


(11a) 


-(n)  .  ^^-(n)  ^  Q 

X 


(n>o) 


(11b) 


-  0<")-np<">  =  0 
X 


(n»l) 


(11c) 


0^”^  +  Jw^"^  =  0 


(n^l) 


(lid) 


-(n)  2-(n)  -(n)  2  r  ,  . 

u  -  ike  V  +  nw  =  —  h  (n>o) 
X  7r  X 


(lie) 


At  x*o  we  want  u,  v,  w,  0  and  p  to  be  continuous.  Writing  fj  =  f(o+)-f(0-) 


we  thus  get: 


-(n)  =  0 

w 


1  w  X  J  “  ° 


(12a) 


(12b) 


using  (11c).  Differentiating  (11c)  and  using  the  other  equations  shows  that 


w  I  ”  w  +  n  p 

L  **  J  L  *  J  L  J 


-J 

L  J 


-(n) 
-n  u 
1  X 

L 


h 

W  X  . 
L  J 


similarly 


=-J 


XX  ■  XX  j 


-n  =  -ikne 


2  !;;<”>  +„2  „(»)  .2ai5^ 

!  X  X  X  .  ’ 


Since  h  »'>^e 


—  i  h 

It  i  XX 


-k|xl-k  /4,  we  thus  get 


(12c) 


=  0 

XXX  ! 


How,  using  (10)  we  can  write  down  at  once  a  particular  solution 
(for  X  4  o)  ' 


-kIxI  „  -k^/4 

*^1^1  2nx  Ke _ _ 

P  =  Pr,e  sgnx  =  — — 5  'Y~2  ^ 7 

(K^-n^+J)  -  e  k  (1+J/k'^) 


e-’^l"!  sgnx. 


(I2d) 


P  being  defined  thereby.  Note  that  (12c)  can  then  be  written  as 


I  zM 

'."xx 


07  2  2  2 

2P  K-n  +J  -  e  k  G+J/K  ) 
o 


(12c) 


P  satisfies  (9)  and  tends  to  zero  at  ±»,  but  is  not  continuous  at  x-0. 
The  solutions  of  the  homogeneous  form  of  (9)  are  combinations  of  exponentials 
where  X  Is  a  root  of  the  quadratic 

X^X^n^+J)  -  eV(x2+J)  =  0 

When  is  small,  one  of  the  roots  is  near  0  and  the  other  near  n^-J;  since 
their  product  is  -e^k^J,  they  are  always  real  and  of  opposite  sign. 


The  solution  of  these  equations  Is 

r  2  2  2  2'  2  2 

A  -  -B  >■  P  i  +e  k  J/K  j/(L  +!.  ) 

C  -  2Pjj  [-L^+e\^J/K^^  /(L^+t^) 

D  -  0 

In  principle,  this  gives  the  solution  to  our  problem,  but  because  of  the 

2  2  2 

somewhat  complicated  way  that  L  and  £  depend  on  k  It  does  not  seem  to 
be  possible  to  Invert  the  Fourier  transform  In  terms  of  simple  formulas. 


(16) 

i 

f 
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However,  we  can  now  use  the  small  parameter  e  to  obtain  more  tractable 
approximations.  We  have,  first. 


L^+  ^  .i<n^-J+eV)^  +  4eVj  '  =  ln^-j|+e^k^(n^+j)/|n^-j|+0(e^). 


and  so 


2  1  *  2  ,  2  -i"  '  (n^-J Liy^  ..  4. 

1  =  -r-  n  -J+  n  -J  9  ‘  > 

2  .  j  ^  L  n  -J|  J 


(17a) 


\  [^J-n^+|n^-Jl  j 


In^-J| 


+  0(e") 


(17b) 


Also, 


-k  -k  74  2 

^  2niT  ^Ke  ,  ,  2. 2  1+J/K^  .  4, 

“ - 2 — 2  1  +  e  k  -2 — 2 —  • 

K  -n  +J 


(18) 


K  -n  +J 


Let  us  write  n^  for  the  integral  part  of  J  ,  supposing  (as  has  already  implicitly 
been  done  in  writing  down  (17)  that  J  is  not  a  square  integer.  Then  equations 
(17)  become 


=  eVj/(J-n^)  +  O(e^) 

-  J-n^+  eVn^/(J-n^)  +  O(e^) 


(n$n*) 


or 


-  n^-J  +  eW/(n^-J)  +  O(e^) 
^}  =  eVj/(n^-J)  +  O(e^) 


I 


(n>n*) 


Using  these  and  (18)  in  (16)  we  find,  setting 

_o  -  -k^/4,.„2  2._,. 

P  ■  2nir  ^e  /(K  -n  +J)  , 
o 

.  -  -O  2,  2  K^_n^+J  ,  „  ^  4,  \ 

A  -  -  B  *  P^€  k  J  - 5-5-  +  0  (e  )  I 


2  2  2 

r(jV)^ 


2^2/  ^  tv2  t2v  / 

C  -2P“(H  >  *  0(tS 

°  (jV)^(K‘‘-n  +J) 


(n<i\  ) 


(19) 


(20) 


(21) 
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2,  2.  4  2  2,  , 

»  «  „o  , .  e  k  (n  -JK  -J  )  ,  4^ 

A  »  -B  =  P  1+  — j  5 — X— +  0(e  ) 

°\  (n  -J)  (K  -n  ^  J) 


C  =  -2P°  e^k^J  .  -r~"-  ■  +0(e^) 

°  k2(  n^-J)^ 


>  (n>n^ ) 


In  order  to  obtain  a  representation  of  the  solution  which  is  uniformly 
valid  over  as  large  a  range  of  x  as  possible  one  should  not,  however,  simply 

-nTxl 

expand  the  functions  e  '  '  and  cosLx  in  powers  of  e.  It  is  better  to 

follow  the  general  idea  of  multiple  scaling.  Introducing  new  variables 

C  “  ex  and  X  =*  e^x  (and  possibly  =  e^x  etc.,  though  we  shall  not  use  these) 

so  as  to  obtain  a  representation  of  the  form  w  =  w^(x,C,X,Y, +  ... 

in  which  the  second  term  is  in  fact  of  order  compared  with  the  first  even 

_2 

when  x  is  as  large  as  e  .  Thus  in  using  the  above  approximations  in  (15) 
we  write  in  the  "wave  case",  n<n*: 


pO 


,  '  -K  X 

+  . , .  e  '  ' sgnx 


°  K'^(J-n  )^ 


4  2  2 

„o„.  ^  ,.2,2  n  -JK  -J 

-P  (1+sgnx)  1+e  k  - x-x — x — x — 

°  (J-n  )  (K^-n  +J) 


1,2  2 

+  "  X 

^  2  3/2  ^ 


"7,  2 

cos  i  fJ-n  X 


(nin^) 


and  in  the  "evanescent  case",  o>n. , 


°  I  K^-n^+J  ,  /  ^ 

„  r  4  2  -2  1  _h/ 

_o  I,.  Z,  2  n  -JK  -J  .  I  \/ti 

-  sgnxP  l+€  k  — 5 - 5 — 5 — 5 -  +. . .  e  I  V 

°  -  (n  -J)  (K  -n  +J)  ' 

2  2 

-P°  (l+8gnx)eVj  coskC(J/(n^-J))’* 

®  K^(n  -J)^ 


... 

-n  +J)  ' 


2(n2-J)3/2 


(n>n*) 
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The  other  components  of  the  flow,  for  n^l,  can  also  be  obtained  by  using 
equations  (11);  the  results  are: 


(26c)  shows  that  for  continuity  of  at  x*o  we  have  u°  =  -4it  ^ 

since  (27)  gives  u°  '  =  -4Tr  ^^e  ^^/(K^-k^e^)  -  2elk|E  we  get 


2n-Sce-^ 


v2  v2  2 

K  -k  e. 


e  -k  /4 
—  e 


2 

.  -K|,|  clkl  -clkllxl 


v2  k2  2 

K  -k  e 


^2  2 
K  -k  £ 


.a.-Ve-'*''*  I  -K 


I  -K  X  -e  k  X 
1  e  '  ' -e  I  ' I  '  j  sgnx 

\  i 


(28a) 


(28b) 


pV-'  =  .  (28c) 

For  small  e,  the  factor  (K^-k^e^)  ^  in  (28)  may  be  replaced  by  K  ^(1+k^e^/K^ . . . ) 

in  the  same  manner  as  in  obtaining  (23)  and  (24).  It  should  perhaps  be  remarked 

2  2 

that  since  these  expansions  in  powers  of  e  are  not  uniform  in  k  -  in  fact 

2  2  2 

they  should  only  be  used  for  e  k  /K  small  -  we  should  not,  strictly  speaking, 

use  them  in  evaluating  the  inverse  Fourier  transforms.  However,  because  of 

-k2/4 

the  factor  e  we  may  cut  off  the  inverse  Fourier  integral  at  a  finite  value 

of  k,  for  Instance  k  =  Ke  ,  thereby  Incurring  only  a  transcendatally  small 

error.  In  the  truncated  integral,  the  expansions  are  legitimate,  and  the 

Integrals  of  its  terms  can  then  be  replaced  by  infinite  integrals  again  with 

only  a  transcendantally  small  error.  Thus  in  fact,  one  need  not  worry  about 

2 

the  non-uniformity  in  k  of  the  expansions. 


It  is  now  relatively  easy  to  invert  the  Fourier  transforms,  using  the 
following: 


2n  J 


(29a) 


15 


1  "  -lkY-k^/4-a|k|  ..  Tr“'^  (3-iy)\  ^  ^  (a+lY)^„  ^ 

e  '  'dk  -  e  Erfc(a-iY)  +  e  Erfc(a+1Y) 

Z  IT  ^  Z  ! 

—00  I 


(29b) 


1  -ikY-k^/4  2  ,  ,  2.-H:  -Y^/(l+16a^)  . 

'  e  cos(ak  +b)dk  =  ir  ^(l+16a  )  e  '' 

-oo  r  2  " 

,  ^  -1,  4aY 

.  cos  *5  tan  4a - +  b 

l+16a^ 


(29c) 


and  “ 

<'  2 

1  -IkY-k  /4 


— ^  —3.^— 

cosak  dk  =  IT  cosh  2aY 


(29d) 


(See,  for  Instance  j^Erdelyl  et  al.  Tables  of  Integral  Transforms,  McGraw 
Hill  1954 j  ;  the  integrals  are  also  not  difficult  to  evaluate  directly.) 
(29b)  can  also  be  written  in  real  form  as 


1  -ikY-k^/4-a|kL.  -Y^  a2 


-h  -Y  a2 

dk  =  ir  e  e  cos2aY  Erfca 

.  /  2 

+  2it  sin2a(Y-n)dn 


In  addition,  for  the  terms,  we  need  similar  Integrals  with  an  extra  factor 
2 

of  k  ;  these  can  be  obtained  from  the  above  by  differentiating  twice  with 
respect  to  Y: 


i  i  e'^^Ve"^'^  =  (Tra3)"^(2-4Y^/a)e"^ 


Ola) 


00  I  2^2 

1.  1  ,2  -lkY-k^/4-a|k|  _  -h  (4a^+2-4Y^)e  ^  e^  co82aY  Erfca 

It  0  QK  “*  iT 


rY  2 


+2n-^^  e^ 


sin2a(Y-Ti)dn  -  4a ir 


2  f  2  r  2  "] 

+  8aYe  ^  je^  sin2aY  Erfca  +  2v  ^  e^  cos2a (Y-n)dn 
L  ^ 

(31b) 


r 


ij  ''^cos(ak^+b)dk  =  ir  "^(l-lba^) 


-*4/ 

2 


-5/4  -Y^/(1+I6a^). 


;-4Y^^~— cos(i4tan"^4a-^2^^^ — -  +b  -8a 


1-H6a2 ,  ^ 


1-H6a^ 


1- 


4Y 
l+16a^ 


aln(*stan^4a 


^  J' 


1-H6a2 


(31c) 


Ilf] 


=  1  (2-4Y^-4a2)co8k2aY 


^a> 


+4aY  sinh2aY 


Old) 


Using  these,  and  setting 


1 

~/(K^-n^-^J) 

IT 

(32a) 

! 

»  X  Qe  ) 

(so  that  F 

o 

i 

8^ 

=  J/lJ-n^l 

(32b) 

v2 

/  ^It  2|-3 
“  4n  jJ-n  1 

(32c) 

i 

* 

_  j  2  ^ 

«  istan"^v|x|-vlxlY^/(l+v2x2)  -f  Vj-n  x 

(32d) 

■jfc.'- 


The  inverse  transforms  of  (23)  and  (24)  give; 

2  2  2  1 

Q'i'f^-’'l’'l8gnx-(l-Ksgnx)(l+v2x2)-V  ^  ^  cos<^  \  (n<n*) 


w^“^-  Q-fe 


(H-v|x|)  sgnx  (n>n*) 
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(33a) 

(33b) 


,2^2. 


(2-4Y^)e''^  "'^1  sgnx-QJ  V  ' 

K  (J-n  ) 


2  2 

K  -n  +J 


,  Y 


'(2-4Y^+46252)e"^^  e^^^^Er£c(BUl)cos20C  +2n~^  8in2$Ula-n)dn 


-4u"SUI  +  SBUIYc"^  e®  ^  Erfc(6|cl)sln2$l5l 

u  r  ^  2  ' 

+  2n”^  e'^  cos2Be(Y-n)dn  '  sgnx 

Jo  j  . 


-Q(l+8gnx) 


4  2  2 

n  -JK  -J 


7  9  2  2 

(J-n  )^(K^-n  +J) 


„  2  2^-5/4  -Y^/(1+v^x2) 


4  i2-4Y^  — cos(t>-2vlxl  i  1-  sln(j»  >  ,  (n<n*) 


l+v^X^ 


^  l+v^x^  i 


(34a) 


„(n)  .  Q  l+J/K  (2-4Y2)e-''  -""'^'sgn  x  - 
2  K^-n^+J 


4  2  2 

.  n  -JK  -J  _ 

— ?  2  2  2  V 

(n  -J)  (K  -n  +J) 


l+vUl/ 


sgnx 


-V(W,gn.)  I  (2-«'-«H^)»sh<26a) 

K^(n^-J)  ' 


+  43E  8inh(2BCY). 


(n>n*) 


(34b) 


(33a)  exhibit8  the  lee  wavee  seen  at  large  positive  x,  the  function  (fi  of 
(32d)  giving  the  phase.  If  X  is  large  enough  that  vX»l,  the  lines  of  constant 

phase  are  approximately  the  hyperbolas 


^2  _  .(J-n^)^[(^X)^-n|  +  (<>-  J)  vX  =  0 
e^v  ' 


which  have  asymptotes  of  slope 


^  -  ±  (J-n^)** 

dx  dx 
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A  picture  of  some  lines  of  constant  phase  computed  from  (32a)  for  the 
case  =  10  is  shown  In  Figure  1.  Of  course,  if  there  is  more 

than  one  lee  wave  mode  (n^>l),  the  asymptotic  solution  will  be  a  sum  of 
such  terms,  and  the  picture  will  be  more  complex. 


The  inverse  transforms  for  the  other  components  of  the  flow  can  be 
obtained  in  the  same  way  from  equations  (25)  and  (28).  We  give  only  the  results 
pertinent  to  u  and  v  up  to  0(e). 


X  K  [ 

y2 

.  (ysin2CY-Ccos25Y)  +  2it  ^e 
Y  2  , 


Y  2  1  '  M 

\  ^  Ycos2C(Y-n)-Csin25(Y-n)!  dn||‘+0(e2)  (35a) 

,<»).  9.1  ' 

n  ]  ^  YJ-n 

i 

K 

.2 


-hx 


.  (Ysin2elelY-6Ulcos2eCY)+2ir  ^e 
1  rY  2 

,  -h  f  n 

-4x  e 

Jo 


Ycos2BC(Y-ri)-8Csln2eC(y-n)  |  dn  -M)(e2) 

J 


1 


(o<n<n*)  (35b) 


(n)  -Y^-k|x1  _ 

“  “  n  1  K  'y-j- 

'•  n  -J 


8lc|cosh26CY-Ysinh28lclY:|>  -KXe^) 


(n>n*)  '  (35c) 
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2  j  2 

e"^  Ye"‘^l^l-e^ErfcUl  acos2CY+Csin2CY) 
IT  K 


-2it 


-h 


Y  _2 


^  ;Ysln2|cl (Y-n)-|c|cos25(Y-n)  dn  '  sgnx+OCe^) 


,<»)  .  22.i»\g„,  ,?I>‘I-62  iia-a  l/'«\ri,=e|?|- 

n  I  (_  K  ^ 


.  rY  2 

•  (Yco823CY+eCsin2B5Y)+2u"^  I  e'^  (Y8in23 1 c|  (Y-n) 

Jo 

)*?  2 

dn  +  (1+8 gnx)  - — 2  (1+V^X2)  5/4  ^ 

-1  J-n 

2  1 
.Ye"^  /<^'’'''^^^\cos(|i-vlX|sln(}))  +0(e2)  ^  ,  (o<n^n*) 

iT^-kIxI  ‘ 

'  'sgnx - 2 —  ® 

n  -J 


2J 

K^+J, 

i 

k2 

K 


.e 


Ycosh26CY-BCslnh26CY)  >+0(e^)  ,  (n>n*) 

\ 


(36a) 


(36b) 


(36c) 


20 


We  now  must  see  how  to  Improve  the  convergence  of  the  sine  series 

2K 

for  w  so  that  It  can  be  used  for  computation.  Since  Qv  -  for  large  n 

nir 

(See  (32a),  and  2/n,  we  see  from  (33b)  that 

(n)  2K  -Y^,  -kIxI  -nlxL 

w  * - e  (e  '  -e  ‘  ')agn  x+ 

o  xn 

+  terms  of  order  n  ^  or  smaller.  (Of  course  the  series  ^*^8ln  nz  is 

ultimately  rapidly  convergent  if  x4o,  but  the  approximation  by  a  moderate 
number  of  terms  becomes  poor  as  x>o.) 


However,  the  series  with  the  asymptotic  form  as  coefficients  can  be 
summed  exactly,  so  by  removing  this  part  we  are  left  with  a  series  which 

00  n 

converges  rapidly  enough  to  be  useful.  The  power  series  E  Z  /n  converges 

to  -log(l-Z)  inside  the  unit  circle  in  the  complex  Z  plane,  the  logarithm 

being  defined  as  0  at  Z=»0  and  with  a  branch  cut  extending  to  the  right 

from  Z-1.  By  Abel's  theorem,  it  also  converges  to  this  same  sum  on  the 

-Ixl+iz 

unit  circle  too,  except  at  Z^l.  Setting  Z«e  this  gives 


?  e-“l*le^”^/n  -  -logd-e" 
1 


...  f  J-log(l-e"'^*cosz-ie-|*'sinz) 

1 


tan 


-1 


n 

■-txl 


T  “  TxT 

l-e  '  ' 


slnz 


cosz 


tan 


slnz 


-cosz 


In  particular,  for  x*o,  we  recover  the  familiar  series 

f  linjLZ.  tan"^  .  tan"W  z/2  =  ^(Tr-z) 

2  n  1-cosz 

Thus  we  have 


(37) 
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Let  u8  set 


2K  -Y^(e'’^l*l-e"”l^n,  (n) 

—  e  .  - +  w 


Then  ve  have 


(o)  2K  -Y^  1  -k|x|  ,  -1  /  slnz 

w  - - e  I  e  '  'ir-z  -  tan  — i — r 


2  le'  '-cosz/J 

*  />  fn) 

+  r  w  sin  nz 
1  ° 

anu  the  series  In  (39)  should  converge  rapidly  enough  for  direct 
numerical  evaluation. 


Investigation  of  the  asymptotic  behaviour  of  the  terms  In  w' 

-3 

(34b)  shows  that  they  decrease  like  n  ,  so  the  series  may  there  be  used 
directly.  (35c)  shows  that  for  large  n 


(n)„  2K  -n|x)-Y^ 
u  "V/  —  e  '  ' 


Thus  we  set 


;(n)  .  _  2K  -Y^-nlxl  (n) 


(nil) 


Since  In  this  case  we  have  a  cosine  series,  we  use 

s  -nixicosnz  -|x|  .-Ix|. 

le  '  ' — —  ■  Re  ^ -log(l-e  '  'cosz-ie  'slnz)^ 


■  log 


and  so 


u  ■  -  —  log  l-2e  '*Uosz+e 

If 


-2|x|  -Y^  .  u^°^.  Sfi(n) 

'  '  e  +  T  +  cosnz 

2  1 


The  terms  v^”^  already  decrease  like  ”2 

n 
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IV  RESULTS  AND  DISCUSSION 


Numerical  evaluation  of  w  and  v  Is  made  for  both  the  single  mountain 
ridge  and  the  mountain-valley  cases,  under  combination  of  the  following 
conditions : 

h  -  1  Km 
o 

e  -  0.2 

J  -  2.0,  7,  13,  and  20 

Based  on  the  formula  for  a  single  mountain  ridge  shown  In  (3),  the  correspond¬ 
ing  formula  for  the  mountain-valley  terrain  case  can  be  written  In  non- 
dlmenslonal  form  as 

-Y^ 
e 

In  which  the  valley  width  or  the  separation  of  ridges  Is  three  non-dlmenslonal 
unit  of  length,  l.e.,  9  Km  In  actual  distance.  For  the  single  mountain  case, 
the  peak  of  the  ridge  Is  placed  at  the  origin.  For  the  mountain-valley  case, 
the  valley  head  Is  placed  at  the  origin,  and  peaks  of  the  ridge  are  at  x«±1.50, 
y»o.  For  the  parameters  used  here,  the  valley  head  has  the  elevation  of  about 
300  m  above  the  horizontal  plain.  The  abscissa  x  and  ordinate  y  In  the 
graphs  concerning  the  vertical  velocity  distributions  are  dimensionless  and 
each  corresponding  unit  Is  equal  to  3  Km  for  x  and  Km  for  z.  The  velocities 
are  also  non-dlmenslonal.  The  vertical  velocity  w  of  unity  corresponds  to 
3.33  msec'  for  the  value  of  the  parameters  used  In  the  present  model.  For 
the  figures  regarding  the  distribution  of  the  horizontal  velocity  component, 

V,  the  ordinate  Y  Is  also  dimensionless  and  Y  ■  1  corresponds  to  y  15  Km. 

The  velocity  v  of  unity  corresponds  to  2/3  msec  The  horlsontal  velocity 
component  normal  to  the  mountain  ridge  Is  not  of  much  Interest  here  and  Is 
thus  not  calculated. 


h 

h 


^-k1x-1.501^  ^-K|x+1.50| 


U  =  10  msec 

K  «  1.25 
H  «  3  Km 
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Figures  2  and  3  show  the  vertical  cross-section  of  the  vertical  •  i 

velocity  distribution  w(x,z),  at  Y*0  over  a  single  mountain,  corresponding 
to  J  *  2.0,  and  the  horizontal  cross-section  of  w(x,Y)  at  t  *  Km  for 
J  »  2.0  respectively.  Both  clearly  show  that  the  flow  is  down  slope  in  the 
lee  of  the  ridge.  Above  the  mountain  ridge  and  In  the  lee  phase  retardation 
takes  place  with  height.  There  is  very  little  retardation,  however,  farther 
away  from  the  mountain,  as  shown  in  Figure  3. 

The  horizontal  distribution  of  w  for  the  case  J  =  2.0  is  shown  in  y 

Figure  3.  The  phase  line  in  the  horizontal  plane  at  the  z  =  1.04  Km  level 
is  concave  toward  the  downwind  direction.  Bow-shaped  or  crescent  shaped 
clouds  will  develop  as  the  manifestation  of  the  vertical  motion  field  where 
there  is  enough  moisture.  Such  clouds  were  reported  by  Abe  (1941).  A  crescent¬ 
shaped  clear  area  that  developed  over  an  overcast  area  with  low  stratus  clouds 
was  recently  reported  by  Edinger  and  Wurtele  (1972).  Figure  4  shows  the  vertical 
cross-section  of  v  at  Y  =  1  (15  Km)  from  the  valley  peak  of  the  ridge. 

-2 

For  Richardson  number  J  =•  10,  the  Brunt— VaisSlM  frequency  is  about  10  , 
which  is  the  approximate  value  for  the  case  of  the  normal  atmosphere.  Thus  a 
smaller  J  means  that  the  atmosphere  is  less  stable.  As  can  be  seen  from 
Figures  5  to  8,  the  phase  retardation  is  more  pronounced  in  the  more  stable 
case  when  J  *  20  (Figure  8)  than  in  the  less  stable  case  (Figure  5).  The  higher 
J  will  allow  high  wave  modes  to  come  into  play.  This  is  clearly  shown  in  these 
graphs.  It  is  of  interest  to  compare  Figures  2  and  5  in  which  the  basic  con¬ 
ditions  are  exactly  the  same,  except  for  the  two  mountains.  The  vertical 
velocity  distributions  for  the  two  mountain  case  is  quite  different  from  the 
single  mountain  case.  The  superposition  of  the  wave  solution  for  each  single 
mountain  at  such  a  valley  width  used  here  leads  to  the  damping.  However,  the 
downslope  motions  still  take  place  in  the  lee  slopes,  and  the  upslope  motion 
in  the  windward  slopes  of  each  mountain.  The  horizontal  wave  lengths  for 
Figures  6  and  8  have  been  shortened  somewhat. 


Next,  let  us  look  at  Figures  12  to  14,  which  depict  the  vertical 
cross-section  of  the  wind  component  parallel  to  the  valley  or  the  ridges 
at  Y  ■  1.  These  results  are  reported  for  the  first  time.  Figures  12  to  14 
show  different  J's,  varying  from  J  =  7  to  20.  A  considerable  Increase  in 
J  will  raise  the  wave  mode.  The  stronger  the  stability  is,  the  more  difficult 
it  Is  for  the  air  to  move  upward.  Consequently  strong  horizontal  motion  will 
go  around  the  mountain.  Above  the  valley,  the  upvalley  flow  (negative  v  in  the 
graphs)  Increases  with  J  and  the  center  of  maximum  v  descends  to  about  the  1.5  Km 
level  at  Y  =  1.  The  wind  component  v  will  change  from  down-valley  Into  up-valley 
at  about  the  600  m  level.  In  the  daytime  near  valley  floors,  wind  usually  blows 
toward  the  up-valley  direction.  Because  the  down-valley  flow  discussed  above 
due  solely  to  the  lee  wave  effect  persists  all  the  time,  it  will  counteract 
the  thermally  induced  up-valley  winds  so  that  the  resultant  wind  becomes  weaker 
than  the  normal  case.  During  the  night,  the  down-valley  wind  attributed  to  the 
lee  wave  effect  will  strengthen  the  thermally  Induced  down-valley  wind.  These 
results  are  in  very  good  qualitative  agreement  with  the  valley  wind  observation 
made  in  Vermont  valleys  which  have  about  the  same  valley  width  besides  a  similar 
wind  vector  at  the  1  Km  level  as  used  in  the  example  here,  by  the  New  York  Univ. 
Group  (Davidson  et  al.,  1958,1963).  They  have  claimed  that  the  case  for  existence 
of  nighttime  down-valley  wind  was  obvious  and  Immediately  convincing  while  the 
case  for  existence  of  its  daytime  up-valley  wind  was  weak. 

Farther  In  the  down-valley  direction  about  10  Km  in  the  lee,  strong  con¬ 
fluent  flow  reaches  the  surface  level  at  Y  *  1.  The  component  of  the  velocity 
V,  in  the  direction  parallel  to  the  valley  or  ridges  in  the  lee  is  closely  related 
to  the  down-slope  lee  wave  motion.  Because  the  flow  in  general  is  down-slope 
and  the  horizontal  velocity  component,  u,  is  maximxim  near  the  upper  portion  of 
the  lee  slope,  the  flow  has  to  force  the  v  component  into  the  down-valley 
direction  from  the  valley  head.  For  a  valley  width  that  is  twice  or  three  times 
the  size  of  the  present  one  (18-25  Km)  used  for  Illustrations,  the  situation  may 
be  Just  opposite,  as  can  be  seen  from  the  distribution  of  v  at  a  position  near  x=7. 
The  maxlimim  up-valley  flow  is  found  less  than  500  m  above  the  valley  floor  and 
the  maximum  down-valley  wind  is  found  at  the  1.7  Km  level.  Some  observations 
reported  opposite  to  that  by  the  NYU  Group  may  have  been  conducted  near  valley 
center  at  a  site  in  a  wider  valley. 
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Some  more  interesting  results  can  be  found  in  Figure  15,  which  shows 
the  distribution  of  the  wind  component  in  the  vertical  Y-z  plane  on  the 
valley  axis  (x  ■  0)  for  J  =  20.  Note  that  near  the  z  =  0.10  level  (*'<  500  m) , 
V  changes  from  down-valley  flow  into  up-valley  flow.  The  flow  becomes 
stronger  going  away  from  the  valley  head  at  Y  =  0,  indicating  that  the  down¬ 
ward  flow  from  the  valley  slope  on  the  upwind  direction  (ridge  at  x  *  1.50) 
turns  toward  the  down-valley  direction  and  enhances  the  thermally  induced 
drainage  wind.  This  Indeed  has  also  offered  an  explanation  as  to  why  the 
maximum  down-valley  winds  are  usually  found  near  the  valley  mouth,  and  the 
height  of  the  maximum  down-valley  wind  in  the  valley  increases  as  the  down- 
valley  distance  increases.  It  appears  that  the  results  obtained  here  are  in 
good  qualitative  agreement  with  observations  and  the  solutions  have  depicted 
the  basic  physics  behind  these  puzzling  and  Interesting  phenomena  reported  by 
Davidson  et  al.,(1958,  1963)  and  studied  analytically  from  another  point  of 
view  by  Tang  (1968) . 
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Figure  3.  The  vertical  velocity  distribution  W(x,Y)  on  a  horizontal 
plane  at  the  level  z  =  £  Km  (same  parameters  and  units  as 

used  in  Figure  2). 
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Figure  5.  The  vertical  cross-section  of  the  vertical  velocity  over 
mountain-valley  terrain  at  Y  =  0  for  J  =  2. 
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The  vertical  cross-section  of  the  vertical  velocity  over 
mountain-valley  terrain  at  Y  =  0  for  J  -  7 . 
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Figure  7.  The  vertical  cross-section  of  the  vertical  velocity  over 
mountain-valley  terrain  at  Y  =  0  at  J  =  13. 
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Figure  9.  The  vertical  velocity  distribution  W(x,Y)  on  a  ho: 

plane  at  the  level  z  =  ^  over  the  mountain-valley 


vertical  velocity  distribution  W(x,Y)  on  a  horizontal 
le  at  the  level  z  =  y  over  the  mountain-valley  terrain 
J  =  20.  (See  units  in  Figure  2.) 
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CONCLUSIONS  AND  RECOMMENDATIONS 


Study  of  the  flow  over  and  around  a  truly  three-dimensional  mountain 
ridge  of  finite  length  and  typical  mountain-valley  terrains  has  been  a 
difficult  task  owing  to  the  complicated  mathematics  involved.  The  present 
study  has  utilized  the  multiple  scaling  approach  to  seek  representative 
solutions  converging  uniformly,  based  on  a  simple  linearized  model.  The 
present  model  is  different  from  previous  studies  in  that  the  solutions  are 
valid  for  large  down-wind  distances  as  well  as  in  the  near  region.  The 
linear  solutions  can  be  superimposed  to  obtain  solutions  for  the  flow  over 
typical  mountain-valley  cases.  Three-dimensional  solutions  have  been 
obtained  which  enable  us  to  understand  better  how  the  three-dimensional 
lee  wave  affects  the  valley  wind  system.  The  solutions  explain  why  a  valley 
with  width  of  about  10  Km  may  produce  a  weak  up-valley  wind  during  the  day  and 
strong  down-valley  flow  during  the  night.  They  have  also  provided  an  explanation  of 
the  increasing  of  the  Intensity  and  the  height  of  the  maximum  down-valley  wind 
during  the  night  as  the  down-valley  distance  increases.  The  results  obtained 
here  are  in  agreement  with  observation  available  and  other  findings.  The  size 
of  the  valley  width  and  the  height  of  the  ridges  lines  are  important  factors. 

In  this  study  we  have  developed  a  model  based  on  a  relatively  simple  assump¬ 
tion  in  which  a  constant  basic  wind  field  is  used  to  avoid  complication  in 
addition  to  the  three-dimensionality  in  the  beginning.  It  is  expected  a 
multiple-layer  analytical  model  with  vertical  wind  shear  and  different  thermal 
stratifications  can  be  constructed  and  can  extend  into  the  higher  atmosphere  to 
include  the  stratosphere.  The  thermal  boundary  layer  interaction  problem  can 
be  attacked  in  a  more  or  less  similar  fashion.  With  the  knowledge  of  truly 
three-dimensional  flow  over  and  around  mountain-valley  terrain  and,  in  addition, 
the  further  analytical  study  to  be  done,  we  can  understand  much  better  the 
Interaction  between  the  stratified  meso-scale  flow  forced  by  the  topography 
and  the  thermally  Induced  boundary  layer  flow.  At  the  same  time,  an  investi¬ 
gation  with  an  Improved  numerical  model  based  upon  the  work  we  have  developed 
earlier  can  be  revitalized  and  modified  to  put  it  into  daily  meso-scale 
prediction. 
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ABSTRACT 


Based  upon  the  contemporary  theoretical  concept  and  the 
incorporation  of  a  model  recently  developed  by  Kazakov  and 
Larriyev  (1978),  a  parameterization  of  the  atmospheric  surface 
layer  interacting  with  active  soil  layer  is  made  by  taking  into 
account  of  precipitation  and  soil  moisture  contributions.  A 
method  for  predicting  surface  layer  meteorological  variables 
is  developed  and  can  be  realized  as  a  subroutine  in  numerical 
prediction  even  over  non-homogeneous  terrain.  The  development 
of  the  proper  heat  and  moisture  fluxes  boundary  conditions  based 
on  this  study  will  facilitate  the  solving  of  meso-scale  and  larger 
scale  circulation  problems. 


I.  INTRODUCTION 


In  order  to  study  meso-scale  atmospheric  dynamics  and 
boundary  layer  problems,  the  equations  governing  the  atmospheric 
surface  layer  and  the  remaining  part  of  the  atmosphere,  and  some¬ 
times  an  equation  of  the  active  soil  layer,  are  usually  solved 
jointly.  Previous  works  that  took  such  an  approach  are  by,  for 
example,  Estoque  (1961)  and  Physick  (1976)  for  meso-scale  dynamics, 
and  Estoque  (1963)  and  Krishna  (1968)  for  atmospheric  boundary 
layer  studies.  Because  of  the  complexity  of  the  problems,  it  has 
been  difficult  to  obtain  the  optimum  realization  of  the  surface 
layer  model  as  well  as  of  the  main  problem,  i.e.,  of  the  model  for 
the  remaining  part  of  the  atmosphere. 

Gutman  (1954)  devoted  his  work  to  soil  temperature  calculation 
where  a  more  convenient  approach  was  applied,  in  which  the  surface 
layer  equations  are  transformed.  The  effects  of  the  surface  layer 
on  the  soil  layer  are  formulated,  although  the  appearance  of  the 
boundary  conditions  tend  to  become  more  complicated.  In  the  present 
paper,  the  advantages  of  this  approach  over  others  is  that  the 
surface  layer  model  and  the  main  problem  for  the  layer  above  can 
be  attacked  almost  independently.  Kazakov  and  Lazriyev  (1978)  have 
recently  taken  such  an  approach  to  parameterize  the  surface  layer 
in  its  interaction  with  the  active  soil  layer  on  the  basis  of 
current  theoretical  concepts  and  observational  data.  In  general¬ 
izing  this  later  model  we  have  taken  into  consideration  the  soil 
moisture  and  precipitation  contributions  based  upon  the  development 
by  Deardor ff  (1978). 


II.  THE  CASE  WHEN  THE  SURFACE  TEMPERATURE  AND  MOISTURE  ARE  KNOWN 


It  can  be  proved  that  the  constant  flux  concept  remains  valid 
for  the  surface  layer  over  topographically  inhomogeneous  terrain 
under  nonstationary  conditions  as  long  as  the  following  assumptions 
can  be  met : 

1.  Characteristic  time-scale  of  the  processes  are  greater 
than  a  few  score  minutes,  and 

2.  Characteristic  horizontal  length  scales  of  the  processes 
are  in  the  order  of  ten  times  the  thickness  of  the  surface 
layer  (i.e.,  '^500  m)  or  higher. 

Under  such  assumptions  the  following  relationships 
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iH 

3Z 


C^Vu, 


K  = 
u  3z 


Cj^W 
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=  C„V(T-T  ), 
eaz  H  g 


K  ^ 
e  3z 


V(q-qg) 


(2) 


Is  =  s'  ;  'i,  =  SS) 

(all  notations  are  explained  at  the  end  of  this  paper) 
can  be  used  as  boundary  conditions  for  meso-meteorological  and 
planetary  boundary  problems  at  z  =  h,  including  the  problem  of 
mountain-valley  circulation  when  the  terrain  slope  is  not  exceed¬ 
ingly  large.  It  is  understood  that  T^,  qg  and  h  in  this  section 
are  given  functions  of  x,  y,  and  t. 


In  the  paper  by  Kazakov  and  Lazriyev  (19” 8)  a  simple  method 
to  calculate  Cjj  and  C^^  as  functions  of  Rig  and  H, 
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was  suggested.  Although  are  not  known  in  advance,  one  can  i 

either  determine  ,  Cj^  and  by  iterations  by  taking  Vj^ ,  and  i 

f 

C„  from  the  previous  time  steps  in  the  actual  computation.  If  one 
n 

wants  to  obtain  the  continuous  fields  of  turbulence  at  z  =  h,  one 

can  use  the  formulae,  ‘ 


and 


i  JLi' 

3  z  I  h 
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^  d  in 


dC 
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(4) 


(i  =  u, 9 )  , 

to  determine  the  coefficients  of  eddy  viscosity  and  diffusivity 
and  their  derivatives  as  shown  in  Eq.  (4).  Empirical  functions 
iti£  and  fj^  are  presented  in  Table  1.  The  calculation  of  the  distri¬ 
bution  of  meteorological  elements  u,  v,  t  and  q  inside  the  surface 
layer  can  be  obtained  from  the  following  equations: 


u  =  u 
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V  =  V 


f 

U  o 


(5) 


T  -  T„  =  (T,  -  T„)  -fAr-r — ^  , 


(6) 


where  is  a  function  of  both  Rig  and  H,  and  can  be  calculated 
by  formulae  as  indicated  in  the  work  by  Kazakov  and  Lazriyev  (1978). 
The  heat  and  moisture  fluxes  at  level  h  can  then  be  simply  expressed, 
respectively,  as 

H  =  pc  C„V,  (T  -  T,  ) 
p  H  h  g  h 

and 

E  -  qj,) 


(8) 


(9) 
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III.  THE  FORMULATION  OF  THE  EQUATIONS  FOR  THE  CASE  WHERE  THE 
INITIAL  SURFACE  TEMPERATURE  AND  MOISTURE  ARE  NOT  KNOWN 


For  determining  the  surface  temperature,  T  ,  and  moisture 

5 

q  ,  the  problem  has  to  contain  the  description  of  heat  and  moisture 
S 

exchange  in  the  active  layer  of  soil,  and  the  interaction  between 
the  soil  and  the  air.  It  is  essential  that  these  equations  are 
simple  enough  for  constructing  the  algorithm  which  must  be  conven¬ 
ient  to  apply  to  the  complicated  problems  of  atmospheric  dynamics. 
First  let  us  begin  the  description  of  the  model  with  the  equation 
of  heat  conduction  in  the  soil. 


3  _  3^9 

T-qr  8  -  V  - 5- 


(z  =  o) 


It  is  well  known  that  the  thermal  conductivity  of  the  soil  depends 
mostly  upon  the  soil  moisture,  and  to  a  lesser  extent  on  temperature, 
the  structure  of  the  soil  density.  However,  the  detailed  information 
concerning  the  soil  moisture  changes  with  depth  and  also  other 
aforementioned  properties  are  not  available.  Thus  it  would  be 
appropriate  to  assume  that,  for  a  first  approximation,  the  thermal 
diffusivity,  v,  depends  on  moisture  only  at  the  ground  surface 
level  so  that  we  can  apply  the  empirical  relationship  to  calculate 
the  soil  diffusivity  as  shown  below: 


0.001  +  0.0004  W^ 
0.27  +  Wg 


(Cm^sec”^) 


which  is  adopted  from  Benoit  (1976)  and  Sellers  (1975). 

In  order  to  solve  Equation  (10)  we  regard  it  as  a  two  point 
boundary  value  problem.  By  assuming  both  the  absence  of  the  heat 
flow  from  the  deeper  layers  of  soil  below  and  also  a  constant  soil 


6 


temperature  at  a  great  depth,  we  may  write  the  lower  boundary 
condition  as 


0|  -  constant  (12) 

for  a  given  location  based  on  available  climatological  data.  If 
we  assume  tliat  the  so-called  viscous-sublayer  is  absent,  it  would 
be  quite  natural  to  acquire  continuity  at  the  temperature  field 
at  the  soil-air  interface,  i.e., 
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z=o 


T 
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(13) 


It  is  noteworthy  that  the  inclusion  of  the  viscous  sublayer  does 
not  actually  complicate  matters  technically  (see  Zilitinkevich ,  1970). 
However,  no  report  in  open  literature  has  indicated  that  the  in¬ 
clusion  of  the  sublayer  in  atmospheric  dynamics  models  has  provided 
better  results. 


As  for  the  second  boundary  condition  at  z  =  o,  we  make  use  of 
the  heat  energy  balance  equation,  as  usual  at  the  ground  surface. 
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The  last  two  relationships  are  suggested  by  Staley  and  Jurica 
(1972)  and  Idso  et  al  (1975).  For  the  soil  moisture  equations 
we  will  adapt  the  force  restore  method  by  writing  the  rate  of 
change  of  soil  moisture  in  non-dimensional  form  as  shown  in  the 
following : 
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^  P  d  '  w 
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=  6400  sec 


d^  ■  10  cm 


d^  =  ^>0  cm 


It  is  seen  from  Eq.  (18)  and  (19)  that  the  rate  of  change  of 
moisture  content  is  directly  proportional  to  t)ie  difference 


between  the  rate  ol  precipitation  and  that  of  evaporation.  The 
inequalities  in  the  brackets  indicate  that  Equations  (18)  and  (19) 


are  valid  only  when  the  values  of  ip  and  41  are  between  o  and  1. 

S 

Beyond  these  limits  the  rate  change  of  soil  moisture  can  be  nega¬ 
tive  after  reaching  the  maximum  value  of  W  that  soil  moisture 

“  max 

cannot  continuously  increase  and  the  extra  moisture  will  be  runoff, 


In  order  to  eliminate  some  unknown  variables  we  use  the  follow¬ 
ing  relationships,  based  on  (6)  and  (7): 


T  -  T  =  (T.  -  T  )  6  (22) 

u  g  h  g 


q  -q 

^a  ^g  ^h  ^g 
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(24) 


From  an  empirical  relationship  suggested  by  Deardorff  (1973) 
we  present  the  saturated  specific  humidity  at  ground  surface  in 
the  following  form: 


qg  =  “Qg 


\  '  V 


(25) 

(26) 
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where 


f  w  /w,  ( w  <  w,  ) 

a'  =  /  S  ,  S 
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(27) 


Q(T)  = 


i:.!!  7.63 

P  I- 


(28) 


Q(T  )  ,  y)  =  lOoU  mb  =  const 


Attention  should  be  paid  to  the  inequalities  in  (26)  and  (27). 

In  case  q  in  (25)  becomes  larger  than  Q  we  should  put  q  equal 
S  S  S 

to  Q  .  Thus  the  possibility  of  condensation  of  moisture  at  grouiid  S'urface 

is  taken  into  account.  We  come  now  to  the  inequality  in  (27), 

which  means  that  after  the  concentration  of  soil  moisture  at  ground 

reaches  a  certain  critical  value  W,  ,  less  than  the  W  ,  the  further 

K  max 

increase  of  soil  moisture  will  not  change  its  influence  on  the  air 
in  the  surface  layer.  The  pressure  in  the  formula  for  saturated 


moisture  and  humidity  in  (28)  is  assumed  constant,  quite  guarantee¬ 
ing  the  required  accuracy  because  this  formula  is  used  in  this 
paper  only  for  one  level,  namely,  the  ground  surface,  where  rela¬ 
tive  percentage  changes  of  the  pressure  are  small. 

Next  vie  shall  discuss  the  question  of  the  formulation  of  the 
initial  conditions  for  (10),  (18)  and  (19).  In  principle  it  should 
be  possible  to  get  experimental  data  about  values  of  soil  tempera¬ 
ture  and  moisture  at  different  depths  of  the  active  soil  layer  for 
a  certain  instant  whicli  we  may  consider  to  be  the  initial  value. 
However,  at  the  protn.nt  such  data  is  practically  unavaiJable  frje 
routine  calculations.  Tlierefore,  we  shall  assume  that  after'  cal¬ 
culations  have  been  made  for  quite  a  long  interval  of  time,  the 


■ 


> 

f 

i 

I 

initial  conditions  becomes  less  important,  because  the  amount  of  • 

moisture  in  the  soil  is  changing  continuously.  The  moisture  change  ! 

f 

will  soon  begin  to  vary  quickly  as  time  goes  on,  depending  upon  the 
future  development  of  precipitation  and  some  external  factors.  At 

I 

the  same  time  the  influence  of  the  initial  condition  of  the  moisture 
gradually  fades  away.  In  view  of  the  minor  role  of  the  initial 
condition,  the  appropriate  initial  condition  under  the  circumstances 
may  simply  be  assumed  as 


W  =  W„  =  o  when  t  =  o 
g  2 


(29) 


i  I. 


Considering  (10),  a  parabolic  equation,  we  solve  it  with  the  so- 
called  differential  equation  with  no  initial  condition  (See  Tikhonov, 
1972).  Following  such  an  approach,  we  can  solve  (10),  (18)  and  (19). 
It  may  be  worthwhile  to  point  out  again  here  that  unless  computations 
have  been  run  for  a  few  hours  the  results  for  the  beginning  time 
intervals  have  no  true  physical  meaning  and  will  not  agree  with  the 
observations  in  the  first  few  hours. 


IV, 


THE  DERIVATION  OF  THE  COMPUTATIONAL  FORMULAS. 


It  is  important  to  construct  an  algorithm  that  will  be  simple 

enough  for  practical  application.  Our  problem  now  is  to  develop 

some  mathematical  expressions  similar  to  (1)  and  (2),  containing 

parameters  which  are  either  known  or  can  be  computed  so  that  these 

expressions  can  be  used  as  boundary  conditions  in  solving  the 

governing  differential  equations  (or  the  atmospheric  dynamics  above 

the  ground.  First,  we  want  to  obtain  the  expression  for  heat  flux 

from  the  soil,  G,  with  known  surface  temperature,  T  .  The  main 

g 


11 


difficulty  is  that  v  depends  upon  soil  moisture  which  changes 
with  time.  However,  the  indicated  difficulty  can  be  overcome 
if  one  can  generalize  the  known  expression  for  the  flux 


1  f 

-7^  J 


dO,  (f) 

C=o 


dt 


which  corresponds  to  the  exact  solution  for  problem  (10)  with 
boundary  conditions  (12)  and  (13)  when  v  is  constant.  When  v 
is  a  given  function  of  time,  t,  (10)  may  be  written  in  the  simpler 
expressions  below,  by  absorbing  v  into  ^  through  a  simple  transfor¬ 
mation. 


80  8  0 


'(t)dt 


The  corresponding  solution  for  (31)  is  then 


ill  -  A. 


ds' 


By  changing  back  C  into  the  original  notation  t,  we  obtained  (33), 
which  is  the  generalized  equation  of  (30)  for  the  case  v  =  v(t). 


H-  =  ^  f  - -  J 


For  convenience,  we  introduce  time  index,  i  ,  and  time  step,  At, 
and  denote 


t^  =  j.At 
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and 


4>(t^  )  =  (j)^  (t  =  T,  q,  V,  etc.) 

Furthermore  we  introduce 

^  j+1  _  Y  1+1  _  ^1  (36) 

For  brevity,  we  shall  drop  the  time  index  hereafter  if  the  time 
step  is  1+1.  We  should  always  calculate  all  variables  for  the 
time  step  1+  1,  as  the  predicted  result  for  the  first  time  step, 
in  terms  of  the  previous  time  steps.  In  order  to  write  (33)  in 
finite  difference  form,  let  us  replace  6^_^(t)  and  v(t)  by  segments 
of  straight  lines  as  shown  schematically  in  Fig.  1.  As  a  result  of 
some  elementary  transformation  we  have 


G  =  X 


90. 

gSz 


=  M  T  +  <5,  A  =  (c-PoV) 
z-o  o  g  \  g  s  s  g 


<5  =  g3  +  N,  t  3  +  N„  t  3  + 

1  g  2  g 


where 


n=m-m3 
1  I.  o  ’  2  2  1 


M 


2  A 


M 


°  ttA  t  V 

2  A  „  1 


/v^  +V  +  "I  V 


M.  = 


2  A 


\itA 


t  Vv^  +  v3  TT  + 


(37) 

(38) 


I 


>  (39) 


The  expression  for  6,  generally  speaking,  is  an  infinite  series, 
but  we  truncated  the  series  at  the  third  term  for  practical  purposes. 
We  have  performed  many  numerical  calculations,  under  the  assumption 


13 


that  V  =  constant,  which  have  shown  that  the  first  three  terms 
can  provide  the  accuracy  as  required  if  the  time  increment  is 
not  more  than  25  -  30  minutes.  Thus  to  use  (37),  it  is  necessary 
to  have  information  of  temperature  at  ground  surface  for  the  time 
steps  j  -  1  and  j  in  order  to  calculate  the  value  at  time  step 
j  +  1. 


As  to  the  computation  of  humidity  or  moisture  at  ground  surface 

for  the  time  step  j  +  1,  it  can  be  obtained  by  solving  (18)  and 

(19)  by  forward  time  difference  method  in  terms  of  values  at  time 

step  j.  Having  G,  one  can  begin  to  transform  the  energy  balance 

equation  (14)  at  ground  surface.  Prior  to  the  transformation,  we 

4  4 

want  first  to  linearize  the  terms  contain  T  and  T  and  Q  ,  by 

g  a  ^g’  ^ 


making  use  of  small  time  step.  At,  in  the  usual  fashion,  i.e.. 


T  =  (T  ^  +  T  )‘^  =  (T  ^)^  +  4(T  ^)^T„  (40) 

g  g  g  g  g  g 


T  ^  =  (T  ^)‘*  +  4(T  ^)^T 


(41) 


and 


(r 


j  dQ(T) 
g  "  dT 


g 

■) 


(42) 


Note  that  because  of  the  linearity  of  relationships  (22)  and 
(36),  the  following  relationship  takes  place. 


(43) 


=  (1  -•  3)  T  +  3t, 
a  g  h 

Furthermore,  if  one  solves  the  system  of  equations  (23)  and  (25) 

for  unknown  variables  q  and  q  ,  we  have  then 

ig  'ia’ 


3  ^h-*-  a(l  -  g) 
3  +  a  -  oB 


(44) 


and 


g 


gQ^  -I-  (1  -  a)  3q^ 
a  +  3  -  a3 


(45) 


.  .  4  4 

Substituting  all  necessary  parameters  H,  E,  G,  T  ,  T  ,  and  R 

g  a 

from  (8),  (9),  (37),  (40),  (41)  and  (16),  respectively,  in  addition 
to  the  estimated  downward  short  wave  radiative  flux  near  the  surface, 
S,  into  the  energy  balance  equation  (14),  and  solving  for  the  tempera¬ 
ture  increment  x  ,  we  have  finally,  after  rearrangement, 
g 


g 


(46) 


where 


LpChVha 


g  +  3  -a3 


and 


A„  =  (1  -  g„)  S  +  e  a-  (T  ^  )‘^  +  6 
T  g  g  R  a 


*  ‘■“pVh  >  -  ''q 


e  a  (T  ^  )‘^  + 
g  g 


(47) 


(48) 


-  P) 
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With  (46),  H,  at  level  A  on  the  top  of  the  surface  layer  at  time 
step  j  +1,  can  be  obtained  from  (8)  as 


H 


H 


H  q 

q  h 


+ 


where 


H 

H 

H 


T 

q 

T 


pc  C„V  A  /A 
p  H  h  T 

pc  C  V,  A  /A 
P  H  h  q 

(pCpCj^V^A^/A)  + 


(49) 


>  (50) 


A  similar  expression  can  be  obtained  for  E  through  a  somewhat  lengthy 
manipulation,  which  is  shown  in  the  appendix,  to  yield 


^  ^  '  h  *  Ih  *  ■ 

where 

E  =  fbA  /A 
E„  =  f(bA  /A-1) 

q  q  ’ 

E^  =  (fbA^/A)  +  Q'^  , 


f 

b 


( _ 2__ 

a+B-a8 


) 


3.80  .63 

P 


4.29  X  10^ 

jol.S6,10^/Tg 


(51) 


(52) 


Except  for  the  coefficients,  variables  in  Equations  (49)  and  (51)  are 
for  the  time  step  j+1  at  level  h  and  are  to  be  determined.  These 
two  equations  will  serve  as  boundary  conditions  for  the  prognostic 
equations  of  diffusion  for  temperature  and  moisture  in  the  atmos¬ 
pheric  layer  above.  This  set  of  equations  can  be  solved  iteratively 
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in  addition  to  the  momentum  equations  and  their  associated 
boundary  conditions,  until  the  solutions  converge  to  a  desired 
accuracy.  When  and  are  determined,  the  values  of  q  and  T 
at  the  surface  and  the  intermediate  levels  in  the  surface  layer 
can  also  be  easily  determined.  The  proposed  more  explicit  step 
by  step  procedure  will  definitely  facilitate  the  numerical  pre¬ 
diction  of  meteorological  variables  in  meso-meteorology . 
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V.  DISCUSSION 


The  horizontal  wind  components,  potential  temperature,  and 
specific  humidity  in  the  surface  layer  can  be  simply  written  in 
terms  of  a  set  of  appropriate  functions  and  corresponding  specified 
boundary  values,  when  the  time  and  length  scales  in  consideration 
are  not  very  small.  Their  fluxes  at  the  upper  boundary  can  also  be 
expressed  in  terms  of  the  appropriate  boundary  values  and  transfer 
coefficients  which  depend  upon  the  Richardson  number  and  the  depth 
of  the  surface  layer.  Thus  this  formulation  and  parameterization 
is  convenient  for  studying  meso-scale  boundary  layer  problems  where 
the  surface  layer  effect  is  important.  The  model  can  be  used  to 
determine  the  surface  temperature  and  specific  humidity  at  the  sur¬ 
face  by  including  the  precipitation  and  soil  moisture  when  the  values 
at  the  top  of  the  surface  layer  can  be  obtained  by  solving  the 
governing  equations  in  the  layer  above.  Prediction  formulae  for 
heat  flux  are  also  developed  by  using  finite  time  difference  formula¬ 
tion. 

This  model  and  parameterization  can  be  applicable  to  planetary 
boundary  problems  as  well  as  to  meso-meteorological  problems  over 
mountain  terrain  (Tang  and  Peng,  1977).  Since  mountain-valley 
circulation  results  from  the  interaction  between  the  air  in  the 
atmospheric  surface  layer  over  the  sloping  mountain-valley  surface 
and  that  of  higher  layers  and  the  ambient  atmosphere,  the  determin¬ 
ation  of  the  meteorological  variables  in  the  surface  layer  is 
essential.  The  formulation  and  parameterization  developed  here  will 
facilitate  the  numerical  experiments  of  three-dimensional,  non¬ 
stationary,  meso-scale  mountain-valley  circulation  and  interaction 
with  the  layer  scale  flow. 
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Note  that  the  earth's  surface,  in  particular  mountain  or 
valley  surfaces,  is  usually  covered  by  forest  and  vegetation. 

Based  on  Deardorff  (1978)  we  made  estimations  which  indicate  that 
the  vegetation  layer  greatly  influences  the  meteorological  processes 
in  the  boundary  layer.  Therefore  it  is  our  intention  to  consider 
the  above  concept  and  take  this  more  sophisticated  approach  to 
improve  the  present  model  of  three-dimensional  flow  over  mountain- 
valley  terrain  step  by  step  in  the  near  future. 
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NOTATIONS 


Arabic 


a  anemometer  level  over  the  earth  surface 
Cg  specific  heat  of  soil 

Cp  specific  heat  of  air  at  constant  pressure 
Cjj  drag  coefficient 

heat  and  moisture  transfer  coefficient 
soil  depth  influenced  by  the  diurnal  soil  moisture  cycle 
soil  depth  influenced  by  seasonal  moisture  variations 
evaporation  rate 


^2 


(  )  the  subscript  g  refers  to  the  variable  in  question  at  ground  surface 

s 

G  soil  heat  flux 

h  surface  layer  heigjit  above  the  ground  surface 

H  sensible  heat  flux  frtan  the  ground  (positive  when  it  is  upward) 

j  tijiB  step  index 

Kp  themel  eddy  diffusivity 

K„  momentum  eddy  diffusivity 

L  latent  heat  of  vaporization 

p  mean  value  of  pressure  at  the  ground 

P  precipitation  rate/nass  per  unit  time  and  area/ specific  humidity 
specific  humidity  at  anemometer  level 
q(T)  saturation  specific  humidity  at  temperature  T 
Rig  bulk  Richardson  number 
R  Downward  longwave  radiative  flux 
S  influx  of  shortwave  radiation  near  ground 
t  time 

At  time  step  interval 
T  temperature 

friction  velocity 

u,v  velocity  ooirponents  in  x-  and  y-  axes  direction,  respectively 
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f 


Arabic 


V 

w 

W2 

x,y 

z 


W, 


K 


W_ 


max 


total  wind  speed 

volumetric  concentration  of  soil  moisture 

vertically  averaged  value  of  volumetric  moisture  concentration  in 
the  majority  part  of  d^  below  which  the  moisture  flux  is  negligible 

horizontal  coordinates 

vertical  coordinate/positive  upward 

roughness  length 

critical  or  saturated  value  of  w 

maximum  value  of  w 


Greek 


g 


g 

e 


w 


degree  of  soil  surface  water  saturation 
fraction 

ground  surface  albedo 
thermal  diffusivity  of  soil 
ground  surface  emissivity 
soil  temperature 
temperature  at  z  =  o 
air  density 
soil  density 
water  density 
Stefan-Boltzman  constant 
cloud  fraction 
buoyancy  parameter 
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Appendix 


Derivation  of  the  working  formula  for  moisture  at  the  top  of 
the  surface  layer,  E. 

From  (8),  (25),  and  (28),  E  may  be  written  as 


^  ct  +  6-oi6  "  ‘’h  i 


A.l 


J 


Because  of  the  non-linear  expression  of  Q  in  T  ,  linearization  was 

o  o 

made  from  (28),  i.e. 

_  3.80  -7.63  I  ’’’S  ~  ■ 


3.80  j^g7.63  jl  - 


242.9®, 


1  ■  Tg 

1 


=  b 


I  10^' 


86x10' 


Tg 


,  ,  3.80  ,n'7-63 

where  =  — —  lo 


Let 


A. 2 


e  = 


1.86x10' 


,  X  =  T, 


10 


1 

X  . 


Expanding  e  into  Taylor’s  series  at  x  =  yields 


A. 3 


1  1_ 

X  X 


1 

X 


1 

X  i 


‘  °  si  V'l^o  "  It  b  '  /x==<„  *  •• 

A. 4 

For  Tg  =  300°K  and  time  interval  of  20-30  minutes,  (x-x^)  will  be 
about  5°C.  in  extreme  cases.  The  higher  order  terms  are  much  smaller 
than  the  first  order  term  and  can  be  neglected.  Then 


1  1_  ^ 

-  e^o=  -e^o  Ax 


A. 5 


23 


!  I 

t  " 


Set  X  =  and  x  =  •  Then 


g 


qJ+I  -  qj  =  b 

g 


1 


1  ^ 


where  b. 


4.29x10^  ^g-1. 86x10  /T 

iii' 


,1^"'  -  ,,bb. 


3  /rn  j 


j  +  l 


-  T  ^ 
g 


A. 6 


and 


3.80 


Writing  in  the  proper 
(46) ,  we  have 


notations  defined  in  this  paper  and  using 


b  b,A 


-  T  +  A  A  E 


bb.  An 


A. 7 


Then  finally  we  may  write,  by  using  (A.l), 

=  ^'h  *  Vh  * 

where 

=  fbb^A^/A 
Eq  =  f(bb;^Aq/A  -  1) 

Eq,  =  abbqA.j,/A  +  • 

and 

OL 

^  h  a+g-oB 


A. 8 


A. 9 


'3 

i 
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